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Abstract 

It is shown that an effective anisotropic spherically symmetric fluid model with heat flow can 
absorb the addition to a perfect fluid of pressure anisotropy, heat flow, bulk and shear viscosity, 
electric field and null fluid. In most cases the induction of effective heat flow can be avoided. There 
is a relationship between anisotropic and charged perfect fluids. 
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I. INTRODUCTION 



Spherically symmetric perfect fluid solutions in general relativity have been studied from 
its very beginning, starting with the interior Schwarzschild solution. Gradually different 
mechanisms in stellar models have been identified that create pressure anisotropy 

a, a. 

Radiating spherical collapse demands the introduction of heat flow [3] or null fluid which 
describe energy dissipation in different approximations [J]. More realistic fluids possess also 
bulk and shear_yiscosity j^], Q], 3]. Charged perfect fluids or dust have been discussed by 



many authors 



It has been shown too that the sum of two perfect fluids, two null fluids 



or a perfect and a null fluid can be represented by effective anisotropic fluid models js|. In 
view of the many existing relations among the fluid models mentioned above it is hard to 
point out some center. 

In this paper it is shown that the anisotropic fluid with heat flow is in some sense the 
most fundamental model and can absorb the addition of viscosity, charge and null fluids. 
Something more, heat flow is not generated in most cases. 

In Sec. II the basic anisotropic fluid model is defined. In Sec. Ill the fluid is supplied with 
bulk and shear viscosity which leads to a new effective anisotropic model. In Sec. IV the 
same is done for the addition of charge and in Sec.V null fluids are accommodated into the 
anisotropic model. Sec. VI summarizes all additions and the effective characteristics of the 
general anisotropic model are given. Several conclusions are drawn. 



II. ANISOTROPIC FLUID MODEL WITH HEAT FLOW 

Einstein's field equations are given by 

8nT Q/3 = G a/3 (1) 

where G a p is the Einstein tensor, T a p is the energy-momentum tensor (EMT) and units are 
used so that c = G = 1. The general spherically symmetric metric is written as 

ds 2 = -A 2 dt 2 + B 2 dr 2 + R 2 (d6 2 + sin 2 d V 2 ) (2) 

where A,B,R are positive functions of t and r only. The spherical coordinates are numbered 
as x° = t,x l = r, x 2 = 8,x 3 = (p. The Einstein tensor involves the Ricci tensor and 
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scalar which are given by the metric and its first and second derivatives |4|. Its non-trivial 
components are Goo, Gqi, Gn, G22 = sin -2 6 G33. 

We are interested in the structure of EMT. It reads for anisotropic fluids with heat flow 

T a p = (/i + p t ) U a Uj3 + p t g a f3 + (Pr - Pt) XaX/3 + QaU/3 + U a q p . (3) 

Here /i is the energy density, p r is the radial pressure, p t is the tangential pressure, u a is the 
four velocity of the fluid (a timelike vector), x a is a um t spacelike vector along the radial 
direction and q a is the heat flux (in the radial direction too). We have 

u a u a = -1, x°Xa = 1, u a Xa = 0, u a q a = 0. (4) 

It is assumed that the coordinates are comoving, hence, the fluid is motionless in them 

u a = A~X, X a = B~ l 5?, q a = qB- 1 5^ (5) 

where q = q(r,t). This gives 

T 00 = fiA 2 , T 01 = -qAB, T n =p r B 2 , T 22 =p t R 2 , (6) 

which should be plugged in the Einstein equations. The anisotropic fluid does not radiate 
when q = and becomes perfect when p r = p t . Thus it accommodates anisotropy of pressure 
and heat flow when one starts with a perfect fluid model. 

Let us see now what happens when other EMTs are added to the basic anisotropic one. 

III. BULK AND SHEAR VISCOSITY 



Bulk viscosity [5j, [6(, [7| adds to the basic EMT the following piece 

1% = -COh a0 (7) 

where ( is a coefficient, is the expansion of the fluid and h a p is the projector on the 
hyperplane orthogonal to u a 

6 = u a . a) h af3 = g af} + u a Uf3. (8) 
Obviously, this means the appearance of effective pressures 

Pr = pf = -ce, (9) 
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which should be added to p r ,pt- They do not change the degree of anisotropy Ap = p r —pt- 
Thus even perfect fluid can absorb the bulk viscosity. The quantities fi, q remain the same. 
No heat flow is generated in particular. 

Shear viscosity is responsible for the piece 

= ^h^h^ 5 (10) 

where a a p is the shearing tensor 

<J a (i = «(a;/3) + «(««/?) ~ (H) 

i] is some coefficient and a a is the acceleration 

a a = Ua-pu 13 . (12) 
The shearing tensor satisfies the conditions 

a aP u p = 0, a af3 g at3 = 0, (13) 

hence, Eq (10) transforms into 



7*, = -2 V a a p. (14) 
Use of Eqs (5,11) gives the non-zero components of the shear J4] 



a n = -B a, a 22 = = --R cr 15 

3 sin v 3 

where 

2 -a 2 = a^a aP . (16) 
One can check that the same components follow when a a p is written as the tensor 

Gap = -TfKp + oXcXn- (17) 

Thus it coincides with the general shear tensor defined by Eq (11) in the spherically sym- 
metric case. It also satisfies relations (13) in any metric. 

Plugging Eq (17) into Eq (14) and comparing it to Eq (3) we find the effectively generated 
pressures 

p s r = -2pf = —rp. (18) 

The degree of anisotropy is changed. There is no generation of energy density or heat flow. 
The scalars 0, a can be expressed through the metric and its first derivatives. 



IV. ELECTROMAGNETIC FIELDS 



The EMT of electromagnetic fields is given by 

T^ M = ± ( - \g aP F^) (19) 

where is the Faraday tensor. One defines a unit timelike vector field n^. An observer 
moving on its direction will measure electric and magnetic field respectively 

F a — Fq^ti^ , H a = —£ afii/F^ . (20) 



These fields are spacelike, E a n a = H a n a = 0. The Faraday tensor decomposes like [1C 

F af3 = e a ^H» - 2E [a np ] . (21) 



Plugging this expression into Eq (19) gives formula (7) from Ref 10j, which becomes after 
some rearrangements 



■■EM _ 1 ( tj2 i 77.2 f ' ^ ' 



T a = ^ [ H + E ) { n * n P + 2 9 ^) ~ 4^ ( E « E f> + H « H f>) + 2 K<* n fi ^ 

where E 2 = E^E^, H 2 = H^H^ and j a is the Poynting vector that measures the energy flow 
in the spacetime 

J a = ^e^R v . (23) 

In order to absorb this EMT by the EMT for anisotropic fluid we choose the direction 
n a = u a and x a — E a / E. The latter is possible because when spherical symmetry is imposed 
H a = and E a has only a radial spatial component. The would be heat flow term in Eq 
(22) disappears and we get 

E 2 

T aB = 2eu a up + eg a(3 - 2ex a Xp, e = ^- ( 24 ) 

<57T 

Thus the addition of electric field induces effective pressures and energy density 

V E = Pf = -V E r = e (25) 

related by simple linear equations of state. Hence, a charged perfect fluid or a charged 
anisotropic fluid may be represented effectively by some neutral anisotropic fluid. There is 
no heat flow induction in this case. 
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V. NULL FLUID 



Null fluid describes dissipation in the free streaming approximation 4| and adds to the 
basic EMT the piece 

where l a is the null vector 



IS = d «h (26) 



T = A~ l 5* + E" 1 ^ = u a + x a , (27) 

satisfying the relations 

1% = 0, Z% = -I. (28) 
Substituting Eq (27) into Eq (26) one finds 

Ta/3 = SU a Up + EXaXp + £ (u a X(3 + UpXa) ■ (29) 

A comparison of this expression with Eq (3), taking into account that q a = qx a , shows that 
the addition of null fluid generates effective energy density, radial pressure and heat flow, 
all of them equal 

fI N =p N = g N = £ (30) 

No tangential pressure is generated. This is the only case where an effective heat flow is 
induced. 

On the other hand, it is known j^, [ll] that a perfect fluid plus a null fluid are equivalent 
to an anisotropic fluid without a heat flow when a rotation of u a , l a is done. This is true 
in any metric. Can we do the same for the sum of anisotropic and null fluid? The Letelier 
method is based on the condition uH^ = 0, imposed on the rotated vectors. Then one of 
them is necessarily timelike, the other spacelike and the EMT of the sum may be written 
as in Eq (3) with q = 0. However, when the basic model is anisotropic, there is already 



a spacelike vector, namely x a ■ Therefore we shall rotate l a into it. Following Ref. [9j we 

P 

(n + p t )u a U/3 + el a l/3 (31) 



change the signature of the metric and extract from T Q/ g + T^l the part 



where now u^u^ = IjX^Xh — ~ 1- The above expression remains invariant under a rotation 
with angle a 

u a = cos a u a + Csina l a , l a = -C _1 sin a u a + cos a F, (32) 



Let us demand that l a is proportional to x a an d thus is spacelike 

l a = ax a , l a l a = -a 2 . (34) 

The relation x >J " u n = yields 

2 sin 2 a 

tan a = Aru^, a = — (35) 

J\ 

Next we find 

sin 2 a 

u^u^ = 1 + sin 2 a, u% = — u%. (36) 

Thus tt a is timelike but is not orthogonal to the spacelike l a . The EMT of the two fluids 
sum reads 

(jj, + p t ) u a up - p t g a/3 + a2 Pt + £ ) 1<JP- ( 37 ) 

Now we normalize u a , l a and then perform another Letelier rotation to make them orthogo- 
nal, uHfj, = 0. This can be done when p r — p t + ea 2 > 0. This holds e.g. when p r > p t . As a 
result u a stays timelike, while l a stays spacelike and the EMT becomes that of anisotropic 
fluid without a heat flow. The expressions for the effective energy density and pressures 

n 

are more complicated than the lengthy formulas in [9( because two rotations have been 
performed, so we omit them here. 

VI. SUMMARY AND CONCLUSIONS 

The results in the previous sections show that viscosity, electric charge and null fluids 
are equivalent to induced energy density and pressures, related by simple linear equations 
of state /j, = np r , n = 0, ±1 and pt = kp r , k = 0, ±1, —1/2. When all such additions are 
combined and absorbed one obtains an effective anisotropic fluid model with 



// = fi + e + e, (38) 

p e r = Pr -(e-- W -e + e, (39) 

Pt = Pt ~ C© + \w + e + e, (40) 

q e = q + e. (41) 



These should be plugged into Eq (6) and hence in the l.h.s. of the Einstein equations (1). 
There are 4 equations for 8 functions (/i,p r ,p t , q, (, r/, e, e). The quantities a, 9 are expressed 
through the metric. One has to impose 4 relations on these functions or set some of them to 
zero in order to obtain a determined system of equations. Several conclusions can be drawn. 

Only the functions of the two modes of dissipation of energy (q, e) have effect upon the 
heat flow. 

Viscosity (£, r/) does not induce effective energy density. 
An important characteristic is the anisotropy factor 

A p e = p e r - p\ = Ap - 2r]a - 2e. (42) 

We see that shear viscosity and charge induce pressure anisotropy. Their absorption by a 
perfect fluid (Ap = 0) makes the latter anisotropic and adds to more sources of anisotropy 
to the usual ones p|. 

Bulk viscosity and null fluid don't induce anisotropy and may be absorbed into an effective 
perfect fluid model with heat flow. 

Charged perfect fluids are related to anisotropic fluids by 

Ap=-2e. (43) 

In addition to the Einstein equations charged fluids have a Maxwell equation but it serves 
as a definition of the charge Q in terms of e. Thus results about anisotropic fluids may 



be carried over to charged perfect fluids. For example, in Ref 



121 ] all static spherically 



symmetric anisotropic fluid solutions are generated by two arbitrary functions, Z and Ap. 



When Ap = their formula reduces to the result for perfect fluids 131 ] . When Eq (43) holds 



instead, we obtain a description of charged perfect fluids in terms of Z and e. 
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